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Abstract. We consider the unoriented two-dimensional Abclian sandpile model on the half-plane with 
open and closed boundary conditions, and relate it to the boundary logarithmic conformal field theory with 
central charge c = —2. Building on previous results, we first perform a complementary lattice analysis of 
the operator effecting the change of boundary condition between open and closed, which confirms that this 
operator is a weight —1/8 boundary primary field, whose fusion agrees with lattice calculations. We then 
consider the operators corresponding to the unit height variable and to a mass insertion at an isolated site 
of the upper half plane and compute their one-point functions in presence of a boundary containing the two 
kinds of boundary conditions. We show that the scaling limit of the mass insertion operator is a weight zero 
logarithmic field. 



PACS numbers: 05.65. +b,11.25.Hf 

1. Introduction 

Since about ten years, logarithmic conformal field theories (LCFT) have been the focus of increasing 
attention. After the first systematic study conducted in T, they have prompted intense work of theoretical 
developments. A key issue is the proper understanding of their representation theory which is considerably 
more complex than in the more usual rational non-logarithmic theories. For these matters, we will refer to 
the review articles [210]; and the references therein. 

One of the reasons for which they have been so much studied is their capability to describe universality 
classes of two-dimensional critical phenomena with unusual behaviours, due to non-local or non-equilibrium 
features. Examples of lattice systems falling in these classes include polymers, percolation, disordered 
systems, spanning trees and sandpile models. 

Sandpile models have been defined by Bak, Tang and Wiesenfeld 0|, and proposed as prime examples of 
open dynamical systems showing generic critical behaviour (so-called self-organized critical systems). The 
specific sandpile model we examine here is the two-dimensional, unoriented Abelian sandpile model (ASM) . 
It is the model originally considered in 0], and remains one of the simplest but most challenging models. 

It is not clear yet whether the ASM can be fully described by an LCFT, but it is perhaps the model where 
this can be most clearly tested, and where the LCFT predictions can be most easily and most completely 
compared with lattice results, making it a sort of Ising model for logarithmic CFTs. In addition the conformal 
fields should have a clean identification in terms of lattice observables. We believe that it is worth pushing 
the correspondence in a concrete and totally explicit model in order to gain intuition for the somewhat 
exotic features of LCFTs. The intricate structure of their Virasoro (or extended algebras) representations 
have direct consequences on virtually every aspect of LCFTs. In particular, the construction of boundary 
states and the interpretation of the possible boundary conditions is an important issue to which a lattice 
point of view can greatly contribute. 

t Chercheur qualifie FNRS 



Pre-logarithmic and logarithmic fields in a sandpile model 



2 



We will specifically focus on certain aspects of the ASM defined on the upper half plane (UHP). We will 
start by extending the analysis of [S] as regards the operator transforming an open boundary condition into 
a closed one, or vice-versa. This operator converges in the scaling limit to an /i = —1/8 pre-logarithmic field 
/Lt, and by looking at its 4-point function, we will conclude that the two channels of its boundary fusion, 

V_i/8X V_i/8-V°P+<, (1) 

must be kept separate, in contradistinction to the bulk fusion. The representations Vo and TZq respectively 
contain those fields which live on an open (Dirichlet) or a closed (Neumann) boundary. In particular, both 
contain the identity field, but with different properties since = 1 while (I) at = 0. 

Our purpose in this paper is to consider various boundary-bulk correlations for ASM observables and 
to compare them with pure LCFT calculations. We focus in this article on local observables only, defering 
examples of non-local ones to future works. An observable to be discussed below is the height one variable, 
already well studied in the literature, and identified in the bulk with a weight 2 primary field. The other is 
the insertion of a unit mass (dissipation) at an isolated site in the bulk of the UHP, which will be shown to 
correspond to a logarithmic field w of dimension 0, the partner field of the identity. By using the boundary 
condition changing field /i, we compute their 1-point functions on the UHP with mixed boundary conditions, 
open and closed, on various stretches of the boundary, and compare them with (numerical) lattice results. 
For the logarithmic field lo, we will use the lattice results to compute its bulk and boundary operator product 
expansions (OPEs) and its expansion in terms of boundary fields. In particular, the usual interpretation of 
the latter suggests that the chiral fusion 7?.o x TZo should contain a channel built on Vq, 

7^o X 7^o = v°p + 7^;5' + . . . , (2) 

and that these two channels should again be kept separate, for the same reason as in The dots stand 
for the other representations which might appear in the abstract chiral fusion, like the representation TZi ^ , 
and which presumably contain the fields living on other types of boundary than open or closed. 

The article is organized as follows. In Section 2, we will give a short description of the Abelian sandpile 
model and its basic features. We will also recall the basics of the c — —2 LCFT that is relevant here, and 
give a summary of the present status of the correspondence ASM/LCFT. 

In Section 3, we investigate the lattice correlations of the boundary condition changing operator fi and 
its fusion Section 4 considers the unit height variable in presence of a boundary which contains both 
open and closed sites, on the lattice and in the LCFT framework. Section 5 and 6 proceed with a similar 
analysis for the insertion of mass at an isolated site, and determine all relevant OPEs and fusions. 

2. Background 

We start by briefiy recalling that model. A more complete account, on this model and on other self-organized 
systems, can be found in the review articles [ZIIHI- 

The model is first defined at finite volume, say on a finite portion L x M of the square lattice. At each 
site i is attached an integer-valued random variable hi which can be thought of as the height of the sandpile 
(the number of sand grains) at i. A configuration C of the sandpile is a set of values {hi}. A site which has 
a height hi in {1,2,3,4} is called stable, and a configuration is stable if all sites are stable. A discrete time 
dynamics on stable configurations is then defined as follows. 

As a first step, one grain of sand is dropped on a random site i of the current configuration Ct, producing 
a new configuration which may not be stable. If C'f is stable, we simply set Ct+i = C[. If C[ is not stable 
(the new hi is equal to 5), it relaxes to a stable configuration Ct+i by letting all unstable sites topple: a site 
with height hi > 5 loses 4 grains of sand, of which each of its neighbours receives 1. Relaxation stops when 
no unstable site remains; the corresponding stable configuration is Ct+i- 

When an unstable site i topples, the updating of all sites will be written as hj — s- hj — Aij for all sites j, 
in terms of a toppling matrix A. If the toppling rule described above applies to all sites, bulk and boundary. 
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the toppling matrix is the discrete Laplacian with open boundary conditions, that is, Ay — 4 for j — i, and 
Aij = —1 if j is a nearest-neighbour site of i. With these topphng rules, the boundary sites are dissipative: 
when a boundary site topples, one (or two in case of a corner site) sand grain leaves the system. This defines 
the sandpile model with open boundary condition on all four boundaries. 

Closed boundary condition is the other natural boundary condition that we may impose. Closed 
boundary sites, when they topple, lose as many sand grains as their number of neighbours, that is, 3 
on a boundary and 2 at a corner, so that no sand leaves the system. As a consequence, the height variable of 
a closed boundary site takes only three values (two at a corner), conventionally chosen to be 1, 2 and 3. The 
rows of the toppling matrix labeled by closed boundary sites are given by Aij = 3 for j = i, and Ay = — 1 
if j is a nearest neighbour of i. 

Thus closed boundary sites and bulk sites are conservative (diagonal entries of A equal the coordination 
numbers so that A^ = 0), whereas open boundary sites are dissipative. One may also make bulk sites 
dissipative j2j by simply increasing the corresponding diagonal entries A^^ from 4 to 4 + , in which case we 
will say that site i has mass U. (Because when there is enough dissipation in the bulk, the sandpile enters 
an off-critical, massive regime, characterized by correlations which decay exponentially, and described by a 
massive field theory ^H], see below.) The height variable hi of a site of mass ti takes the values between 1 
and 4 + t,j. In this sense, an open boundary site is a closed boundary site with mass 1, but one could also 
consider boundary sites with larger masses. 

In all generality, a specific model is completely defined by giving the values of the masses, in terms of 
which the toppling matrix can be written 

{ti + coord, number of i for i = j, 
(3) 
— 1 if i and j are n.n., 

with ti > 0. Stable sites have height variables hi between 1 and A^^, unstable sites have hi > An + 1. In all 
cases, the toppling updating rule hj — s- hj — Aij applies. Then the dynamics described above is well-defined 
provided that not all sites are conservative, i.e. ti > for at least one site. 

Dhar 10 made a detailed analysis of these models. Under very mild assumptions, he showed that 
on a finite lattice A, there is a unique probability measure on the set of stable configurations, that is 
invariant under the dynamics. Moreover, is uniform on its support, formed by the so-called recurrent 
configurations. The number of recurrent configurations, which plays the role of partition function, is equal 
to the determinant of the toppling matrix, Za — det Aa. One is interested in the infinite volume limiting 
measure P*, defined through the thermodynamic limit |A| — > oo of the correlations of the finite volume 
measures. In this article, we will consider A going to the discrete plane I? or half plane Z x Z>. 

The infinite volume limit will of course depend to some extent on the mass values ti. It is however 
widely believed that the continuum limit of these measures are c = — 2 logarithmic conformal field theoretic 
measures, possibly perturbed. 

The c = — 2 LCFT is the simplest and the most studied of all logarithmic theories. The theory which 
is thought to be relevant for the ASMs is the rational theory defined in JTj, equivalently the local theory of 
free symplectic fermions (12) defined by the action 

s = - [ deBe. (4) 



It has a W-algebra built on three dimension 3 fields, which organizes the chiral conformal theory in six 
representations. Four of them, V_i/8, Vo, V3/8 and Vi, are highest weight irreducible representations 
constructed on primary fields, while the other two, TZq and TZi, are reducible but indecomposable, and 
contain respectively Vq and Vi. Of special interest here are the representation V_i/8, with highest weight 

field /z, non-local in the fields 9, 9, and the two representations Vo and TZq. Vq has the identity I as primary 
field, and TZq has two dimension fields as groundstates, the identity I and its logarithmic partner tu ^ ■.90:. 
The corresponding boundary LCFT has been discussed by several authors ^H], and reviewed in |T^[T^ . 
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Various lattice correlations in the ASM have been computed in order to probe the adequacy of the 
description by a c = —2 LCFT. The most natural and simplest observables, but by no means the only 
ones, are the height variables, namely the random variables given by 5{hi — a) for a = 1, 2, ... , A^^. In this 
regard, the height 1 variable, and those for a > 4 as the case may be, is much different from the other three 
[a = 2, 3, 4), and is much easier to handle. 

On the infinite plane K = 1? , the power law (~ r^**) of the correlation of two height 1 variables was 
established in ("IS, , for the case where all bulk masses are set to zero, while its exponential decay was proved 
in when the bulk masses are all equal and non-zero. This was further investigated in ^Hli where the 
scaling limit was directly computed for the mixed correlations of unit height variables and of an other dozen 
cluster variables |, for all bulk masses equal (to zero or not zero). The scaling limit of the unit height 
variable, among others, was identified with a dimension 2 field made up of simple combinations of lo and of 
derivatives of 0, 0, for the Lagrangian perturbed by a mass term m? J ui. These identifications have been 
recently confirmed by the calculation of all multipoint correlators, and extended to other local observables 
jl9|. The unit height variable in presence of an infinite line of massive sites, crossing the whole plane, was 
considered in PU] . 

A number of calculations have been done on the upper half plane A — Z x Z> (or equivalently on an 
infinite strip). The bulk 1-point function of the unit height variable, in presence of an open or a closed 
boundary, has been worked out in |21| . and is consistent with the field identification mentioned above. The 
boundary 2-point correlators for all pairs of height variables have been computed in [221 again for both 
an open and a closed boundary. More recently boundary 3-point correlators were obtained in |2S|) and 
independently in |24| . where the same computations were carried out in the massive regime (all bulk sites 
have equal mass). Reference |23) also identifies the insertion of a unit mass (or dissipation) at an isolated point 
of the boundary. These results are all consistent with an LCFT interpretation, for a specific identification 
of all boundary height variables in terms of 9, 9, both in the massless and in the massive regime. 

All together, these results provide a lattice realization for some of the fields in the representation TZq, 
in terms of height variables. A lattice interpretation for the (chiral) primary field /z of the representation 
V_i/8 was given in [5] as the boundary condition changing operator between an open and a closed boundary 
condition. More on this identification will be given in the next section. 

In what follows we reconsider the height 1 variable and we discuss the insertion of dissipation at an 
isolated site in the bulk of the half plane, when both the open and the closed boundary conditions are 
imposed on different segments of the boundary. This allows us to discuss several different issues at the same 
time: the boundary fusion of the field /i, mixed boundary/bulk correlators and the boundary and bulk fusion 
of the field to. From now on, we use the conventional model, namely all sites are conservative except: (i) 
those on an open boundary which lose one sand grain upon toppling, and (ii) the few bulk sites which receive 
a unit mass. 

3. The boundary condition changing field 

As explained in the previous section, open/Dirichlet and closed/Neumann are natural boundary conditions in 
the sandpile model, and to our knowledge, the only known ones. According to general principles of boundary 

I The results of pertaining to the cluster variables called Sio and Sii do not refer to the clusters pictured in Fig. 1 of 1181 . 
which are not weakly allowed configurations in the sense defined there. Instead the results for 5io and 5ii reported in Table I 
should be divided by 4, and refer in each case to any of the following three clusters 



Also in Eq. (6.2), there is a missing factor P(l)'^, the minus sign in front of M^'/ld must be suppressed, and there is a missing 
"c.c." in the third parenthesis. There are two misprints in (4.8): the denominator of A/^ should be 8 and not 4, and the two 
terms multiplied by ^/^/tt^ should be separated by a + and not a — . P.R. thanks Frank Redig, Shahin Rouhani and Monwhea 
Jeng for discussions about the clusters Sio and Sn, and Monwhea Jeng for pointing out the misprints. 
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conformal field theory [25) . the change from one boundary condition to the other one is implemented by the 
insertion of a boundary condition changing field, which is also the ground state of the cylinder Hilbert space 
with the two boundary conditions on the two edges. That chiral field /i was determined in IS', and shown to 
be a primary field of conformal weight —1/8. We first recall the analysis which established this result, and 
then extend it in order to discuss its boundary fusion. 

Closing sites on an open boundary or opening sites on a closed boundary changes the number of recurrent 
configurations of the sandpile model, and thus the partition function. Let us denote by Zd{I) and Zn{I) 
the partition functions of the sandpile model on the UHP with open resp. closed boundary condition all 
along the real axis, except on the interval I where the sites are closed resp. open. For 7 = 0, the boundary 
is either all open or all closed. 

The effect of opening or closing sites on / can be measured from the fraction by which the number of 
recurrent configurations increases or decreases, i.e. from the ratio of partition functions Z £>{!)/ Z and 
Zn{I)/Zn{^)- These ratios correspond to the expectation values of the closing or the opening of the sites of 
/, and should be given in the scaling limit by the 2-point function of the two fields located at the endpoints 
of /. For what follows, it is worth recalling some details on the way these quantities are actually computed. 

The partition function of a sandpile model is equal to the determinant of its toppling matrix. Because 
opening or closing sites changes by 1 the relevant diagonal entries in that matrix, the ratios of partition 
functions are given by 

Si = ^^^Iri^ = f + = ^'■'<' + ^■''>«' ■ 

In these expressions, Aop and Ad are the usual discrete Laplacians on the UHP with either open or closed 
boundary condition on the real axis, and {Bi)ij = Sij S(i € I) is a defect matrix which is used to insert or 
to remove a unit mass from the sites of /. So Aop — Bj is the Laplacian on the UHP with open boundary 
condition except on the interval I which is closed, whereas Ad + Bj is the Laplacian for a closed boundary 
condition except on I which is open. 

The first ratio is simpler. The determinant det(I — A^p)ij^j has finite rank and is finite since 
the entries of A~p'^ are finite. Therefore, although the two partition functions are infinite, their ratio is a 
finite number. It is not difficult to see that this determinant has the Toeplitz form, and can be evaluated 
asymptotically when its rank is large, by using a generalization of the classical Szego theorem. For an interval 
I = [l,n], one finds, after subtracting a non-universal term related to the change of boundary entropy, a 
ratio Zd{I)/Zd{9) ^ An^^^ for large n, for A = 1.18894. As this should be equal to a 2-point correlator, 
one finds (/Li^'^(0)^^'-^(n)) — An^^'^, where the constant A is a structure constant of the n field. 

The calculation of the ratio of partition functions for the converse situation is in principle similar but 
brings a little though crucial difference. The ratio is again given by a finite rank determinant det(I+A^^)i.jg/, 
but unlike the previous case, it is infinite because the entries of A~j^ are infinite (see the Appendix). As 
the entries of the matrix I + A~j^ all contain the same singular term 2Apj^jjg{0; 0), its determinant has the 

form 2Ap[!|^jjj,(0; 0)/(n) + g{n). The regularized ratio Zn{I)/ Zn{%) is defined to be the function /(n), and 
can be seen as the partition function Zn{I) normalized by Zn{{0}) = 1 + A^^(0; 0) = 3/4 + 2A^^^^^{0; 0), 
the partition function for a closed boundary with a single boundary site open. 

The regularized determinant itself is given by a determinant. Set Mij = (I + A^j^)^^/. We define a 
new matrix M' from M by first subtracting the first row of M from the other rows, and then subtracting 
the first column of the new matrix from the other columns. The resulting matrix M' has by construction 
the same determinant as M, but has finite entries, except M[i which is the only one to contain the infinite 
term 2A-i„jO;0), 

M{i=Mn, Mlj- - A/y-Mn, (j > 1) = Af,i - Mii, > 1) (7) 
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M^. = M,, - Afy - Ma + A/n. (», i > 1) (8) 
Thus the coefficient of 2Apj^,^^(0; 0) in detM is equal to the (1,1) minor of M', 

f{n) = det(I + A;/)rcg = det(A/')»j>i- (9) 
The regularized ratio Zfj{I)IZj^{^') can be computed exactly when the size of / = [1, n] gets large. One 
finds exactly the same result as in the open case, namely that Z (I) / Z (j!)) ~ An^l"^ for large n, with the 
same value of the constant A. As expected, it implies the same 2-point function {^'^'^ {0)ii^'^ (n)) — An^l^ 
as before. 

These results suggest that the field changing a boundary condition from open to closed or vice-versa is 
a chiral primary field of weight —1/8. The way they fuse on the boundary depend on their relative positions. 
The OPE yP (^^[i^ (£) must close on fields that live on a D boundary, i.e. on fields in the representation 
Vo, while the other OPE ^^'■°(0)^^'^(z) should only contain fields that live on an N boundary, and therefore 
taken from the representation TLq jH] . Hence one infers that 

M^'^(z)m^'^(0) = zV4^i^j^(0) + ... (10) 

^^■^(z) A^^^^(O) = [Aat Ijv(O) log z + c^jv(O)] + . . . (11) 

where subscripts have been added to stress the type of boundary the fields live on. The coefficient \m 
specifies the transformation 

ujn{z) + \n log ^ (12) 

under a chiral conformal transformation z — > w(z). 

For these OPEs to be consistent with the 2-point functions given above, one should have, for a choice 
of normalizations, 

(Iz,)op = l, Gl^^^^A. (13) 

(lAr)ci=0, (c^iv)ci = l, G'^^..^.^^- (14) 
These equations emphasize the fact that Id e Vo and Iat G 'R-q are genuinely different fields, and that the 
two conformal blocks in the /i/i OPE must be kept separate. This is in contrast with the fusion in the bulk, 
where the two identities have to be identified, so that the fields in the Vq channel are to be considered as a 
subset of those in the TZq channel. 

Let us also note that in the lattice ASM interpretation, the limit lim^^o ^4P74 (M^'^(0)/i^'^(-z)) = 1 
should correspond to the expectation value of the identity in presence of a closed boundary condition, and 
that this is actually obtained by selecting the lo channel in the OPE. One may see this as the trace left by 
the regularization used in the computation of the ratio of partition functions, or equivalently by the single 
open site left in the otherwise closed boundary. It fits the picture we are going to give in the next sections, 
where the insertion of a unit mass at an isolated closed site corresponds to the insertion of a field uj (see 
also 122) • By supplying the same arguments to a finite (in the continuum) interval [0, z], one could see the 
opening of that interval as the insertion of a massive defect line, 

M^'^(0)Ai^'^(z) - / Axujn{x). (15) 

Instead of closing or opening sites on a segment, one can do it on two segments, 1\ = \z\,Z2\ and 
^1 = [•23,2:4], a situation that was only briefly discussed in [5]. In the scaling limit where all distances 
I Zj — are large with finite ratios, the appropriate ratios of partition functions Z£i(/i, /2)/^_d(0) and 
Zjv(/i, /2)/^Af(0) (regularized) are expected to converge to the appropriate 4-point function 



ZD{h,h 



Zd{^) 



il _^ (;X^'^(Z1)M^-^(22)M^-^(23)M^'^(24)) , (16) 



^ (^^^^(zOM^'^(22)M^-^(23)A^^^^(24)) . (17) 
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The irreducible representation V_i/g being degenerate at level 2, the 4-point function of /i satisfies a 
second order differential equation. Introducing the cross ratio of the four insertion points x ~ f^fffj with 
Zij — Zi — Zj , one finds the general form T 

{^x{z^)^l{z2)^l{z:i)^x{zi)) - {z^2Zzif'^{l-xf/^[aK{x)+l3K{l-x)], (18) 

in terms of the complete elliptic integral 

dt f f + . . . for a; 0+, 

Jo y/l-xsm^t \-^\og^ + ... fora;-l . 

If we assume the ordering zi < 2:2 < 23 < 2:4, the cross ratio x is real and positive in [0, 1]. 

To verify the asymptotic behaviours (|16|l and H17() explicitly, the coefficients a and (3 must be fixed in 
each case. 

In the first case, namely two closed segments in an open boundary, the ratio 
Zpjluh) det [Aop - Bi, - Bi,] 1 
Zn{$) = = ^''[^ ^2°) 

is a determinant of dimension + |/2|. It contains a 'disconnected' piece proportional to the product 
[Z / Z d(^)\[Z oih) / Z d{^)\i plus connected contributions that involve off-diagonal entries A~p{i,j) with 
i in Ii and j in I2, or vice- versa. These entries decay like the inverse distance between i and j, a distance 
bounded below by 12:32!, implying that the off-diagonal blocks go to when this distance goes to infinity. 

As the limit Z32 going to +00 means x going to 0, the previous constraint fixes unambiguously the 
constants to a = 2^ and P — 0, yielding 

942 

(m^'^(zi)m~-^(z2)m''^^(^3)m^^''(^4)) = — (^12^34)^/^ (1 - xy^^ K{x). (22) 

TT 

If the lengths of Ii and I2 are large enough so that Zd{Ii)/Zd{^) and Zd{I2)/Zd{^) can be 
approximated by A|zi2|^/* and A|2;34|^/^, then the intended verification of lfTC|) can be restated quite 
concretely as the statement that the following convergence holds 

Zu{h,l2)Zum ^ det [I - A-p^ i?,, - B.J ^2 
ZD{h)ZD{h) det[I- Aop'B/J det[I- Aop^B/J tt ^ ' ^ ^ ' 

Unlike the determinants for one interval, we have not been able to do an exact asymptotic analysis 
of the determinant for two intervals. The ratio of lattice determinants (all well-defined and finite) have 
been computed numerically and plotted in Figure 1 as colour dots (details on the numerical aspects of the 
calculations are given in the Appendix). We have taken two segments of equal and fixed length n, lying N 
sites apart. For fixed n (taken to be 30 and 50), the distance N = Z32 was varied from 1 to 80, and the 
numerical results plotted as a function of a; = ( „"jy )^, which ranges between and 1 as varies. The 
scaling regime corresponds to large n and large N., and within this setting, is best approached for large 
enough values of A^, that is for values of x not too close to 1. The CFT prediction for the corresponding 
quantity, namely the function on the right side of Ea. l|23|) . is plotted as a solid line. 

The agreement is more than satisfactory in the region closed to the scaling regime, and when x is 
close to 1, it improves with larger values of n. In particular, this supports the view that the correlator 
{n^'^ ^^'^ 11^'^ ^^•^) is regular at a; = 0, and the absence of the conformal block related to Ar(l — a;). 

The OPEs quoted earlier in this section readily follow from this correlator. When Z12 goes to 0, x goes 
to as well, and the expansion yields 

(M^'^(zr)M^'^(z2)M^^^(z3)M^'^(^4)) = A' zli' zli' + . . . , (24) 
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Figure 1. Ratio of the ASM partition functions in Eq.((23J for two segments of closed sites 

and [n -\- N ,2n + N]) in an open boundary, as function of the anharmonic ratio x. The solid curve 

is the CFT prediction, the other two are numerical: n = 30 (orange) and n = 50 (blue). 



as expected. 

If ^23 goes to 0, then x goes to 1 where the same correlator develops a logarithmic singularity, and we 

find 

/ D,N I \ N,D( \ D,N/ \ N,D/ \\ ^ 1/4 



'23 



1/4 



\0gZ23+zl{'^ log 



Zl4 



I6Z13Z34 



(25) 



which implies (Aa? is the value of the off-diagonal entry in the rank 2 Jordan block) 
Xn = • 



(26) 



Let us briefly discuss the second situation, with two open segments in a closed boundary. The relevant 
correlator {jj.^'^ ^^'^ jj.^'^ ^^'^) can be obtained from the previous one by a cyclic permutation of the 
insertion points, e.g. the permutation (1,2,3,4) — > (2,3,4,1) (it can also be obtained from an inversion, 
a method we will use in Section 4). This changes x into 1 — a;, the prefactor (^12234)^^^ (1 — x)i/'' stays 
invariant, and one readily obtains 



2^2 



{^Ji^'^{z,)^Ji^'^ {Z2)^^^^^{Z3)^^^^^ {z,)) = — {z^.z^,)'/^ (i - xf/^ k{i - x). 



(27) 



It makes use of the other conformal block, and as a consequence, is now logarithmically divergent at a; = 0. 
Proceeding as before, one is led to verify the convergence 



ZN{hj2)ZN{%) 



-{l-xY'^K{l-x). 



(28) 



The ratio of regularized determinants has been computed numerically, in the same setting as for the 
other situation, namely two open segments of length equal to 30 (orange in Figure 2a) and to 50 (blue) . We 
have included larger distances N between the two intervals, namely TV running from 1 to 200 sites, so that 
the variable x could assume smaller values. In order to approach even better the small x region, where a 
logarithmic divergence is expected, we have also considered slightly shorter intervals, of length equal to 20, 
but separated by larger distances, up to 410 sites, which allows a minimal value of x equal to 0.0022 (see 
Figure 2b). 
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Figure 2. Ratio of the ASM partition functions l)28|l for two segments of open sites and 
[n + N,2n + TV]) in a closed boundary, as function of the anharmonic ratio x. In both figures, the 
sohd curve is the CFT prediction, given on the right side of Eg. 1281 . The colour dots show the 
numerical results: in Fig. 2a, n = 30 in orange and n = 50 in blue, where, in both cases, is 
varied from 1 up to 200 sites; in Fig. 2b, n = 20 with A'^ running between 50 and 410. 



4. The unit height variable 

In this section we examine a first instance of a correlator including boundary and bulk operators. The lattice 
quantity we consider is the probability Pi{z) that a certain site z in the UHP has a height variable equal 
to 1, given a boundary condition on the real axis mixing both open and closed sites. 

On a finite lattice A, it was shown in 16 that the number of recurrent configurations with — 1 is 
equal to the total number of recurrent configurations of a new sandpile model. The new model is defined 
from the original one by cutting off the bonds between the site z and any three of its four neighbours, and 
by reducing the toppling threshold of z from 4 to 1, and the threshold of the three neighbours from 4 to 
3. Thus the toppling matrix of the new model is equal to AJ^™ = Aa — Bz, where the matrix Bz is zero 
everywhere in A except at z and the three neighbours, where it reads (first label corresponds to the site z) 



Bz = 



( z -1-1 -r 

-1 1 

-1 1 

V-i 



(29) 



Since det Aa is the total number of recurrent configurations in a sandpile model of toppling matrix Aa, 
the probability P\(z) is given by the ratio det AJ(''^/ det Aa. In the infinite volume limit A UHP, one has 

a 4-by-4 determinant. Here A is the Laplacian on the UHP subjected to whatever boundary condition we 
choose on the boundary. 

The boundary conditions we want to consider are the same as in the previous section: an open 
boundary with closed sites on a finite interval / = [zi,Z2], or vice-versa, a closed boundary with open 
sites on / = [zi, Z2]. We define two probability functions relative to these two situations, P°^{zi, Z2] z) and 
Pf'(zi, Z2; z). Their limit when 3z — )■ +00 yields the bulk probability for any given site to have height 1, 
Pi = P°P(zi,Z2;+zoo) = Pfi(zi,Z2;+zoo) = 2(7r - 2)/7r3 (TBI. 
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The lattice calculation of the two probabilities is straightforward. To the two boundary conditions 
correspond the Laplacians Aop — Bj or Ad + i?/, where {Bj)i j = Si j 5{i e /) is the matrix used in the 
previous section, and one obtains 

^oP(, ^, . det[I-A,-Mi?, + i3,)] 1. _ det[I + A-^(i?,-i3,)] 

det[I-AoVi?,] ' det[I + A-/i3,] ' ^^'^ 

One may take I — [—R, R\ centered around the origin, so that the probabilities are the ratio of a rank 
2R + 5 determinant by a rank 2R + 1 determinant. We note that although the entries of contain the 
divergent piece we called 2A~j!^jj^(0; 0) in Section 3, the ratio defining P^^{zi,Z2] z) is well-defined and finite. 
If / is not empty, the two determinants are proportional to 2/S.^^^^^{Q] 0) (the proportionality factor for the 

denominator, a function of zi — 2:2, goes to Az^^ in the scaling limit), and the ratio is finite. If / is empty, 
the probabihty Pf{z) reduces to det [I — Ajj^ Bz], which is finite because the row and column sums of Bz are 
zero, with the consequence that A^"'^ Bz depends on differences of A~j^ entries only, which are well-defined. 

Precisely in the extreme case when the interval / is empty, the two probabilities Pi^(z) and Pi^{z) for 
having a height 1 at site z in front of an all open or an all closed boundary are given by 4-by-4 determinants. 
For z — X + iy, the expansion of the two determinants in powers of y yields |21| 



P°^{z) = P, + ^ + ..., P^'(z)=A-A + ... (32) 



The two functions P^^{—R, R; z) and P{^'(— i?, i?; z) have been computed numerically for various values 
of R and z. The results for P'^^ {—R, R; z) along different lines in the UHP are pictured in Figure 3. 




Figure 3. Values of — — - — ^^■^+'^) — 1 as function of y (left) or as function of x (right), either 
computed numerically (dots) or from CFT computations, Eq. 13911 (solid lines). The curves on the 
left correspond to 2; = and 5 < y < 200, for values of R equal to _R = 10 (light blue), _R = 30 
(blue) and _R = 50 (violet). Those on the right correspond to i? = 30 and y = 15 (light orange), 
1/ = 40 (orange) and y = 70 (red), with -100 < a; < 100. 



If (j){z) is the scaling field to which the random lattice variable 5{hz — 1) — Pi converges in the scaling 
limit, P°^{zi, Z2; z) should be given by the following correlation 



{t,^'^{zi)^i^^^{z,)^{z)) 



^1 (zi, Z2; z - — nrwi — ^ — WW — • ^"^"^i 
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The 3-point correlation function on the UHP can be rewritten as a 4-point chiral correlation on the whole 
plane 

P°P(^l, Z)-Pi = il^'^'^izi) ^^"^^^(^2) </>(z) ^(Z*))chiral,planc • (34) 

Az^2 

In concrete terms, the field is a primary field of weight (1,1), proportional to dOdO + dOdO in the 9, 6 
model [^. The relevant 4-point function is then easy to compute. It satisfies two second order differential 
equations: one coming from /i, the other from <j). By combining them, one obtains a first order differential 
equation, whose solution yields the general form of the chiral 4-point correlator {x = f^fffj) 

1/4 D 

(^(zi) ^l{z2) (f>{z3) (j>{z4)) = a (35) 

The constant a can be fixed by requiring that the expression 

P°P(zi, Z2; z) - Pi = --^ (36) 



4Ay2 



X 



reduces to +Pi/iy^ when zi goes to Z2- 

Clearly 212 going to corresponds to x going to 0. However, Z12 can go to in two ways: either zi 
and Z2 meet at a finite point of the real axis, or they meet at infinity, by going "around" the equator of the 
Riemann sphere. In both cases, x approaches 0, but from different directions. To see this, we set zi = —R 
and Z2 — R, and consider 1 — x for Z3 = z^ = z 

Z14Z23 {R- z){R + z*) 
I — X ^ = n r- (37) 

213Z24 {R-z*){R + z) ^ ' 

This expression makes it manifest that 1 — x has unit modulus and circles around when R ranges from 
to infinity. The expansion of x around R — and R = +00, 

--— + ... forP^O, 

. 1"^! (38) 

^ + ... ior R^ +00, 

R 

shows that x draws out a unit circle centered at 1, travelled anticlockwise from — iO to +iO as R ranges from 
to +00. 

One can write x ~ 1 — e'^ with 6* ~ 0"^ for i? ~ O"*", and 9 ^ 2tt^ when R ~ +00, implying that 
limfl^o Vl — X = +1 and limfl-j+oo \/l — x = — 1. Therefore when expanding the correlator (|36|l around 
X — 0, one takes the positive square root — x when R is close to 0, and the negative root — vl — a; 
when R is close to infinity. 

Bearing this in mind, the expansion appropriate to recover the fully open boundary is around i? = 0, so 
as to shrink the closed portion to nothing. This fixes a = and in turn, yields an explicit expression 

for the probability in the scaling limit 

Pi P2 _ U|2 

P-(-P,P;z)-Pi = --^il^. (39) 

This expression is compared with a lattice calculation of the probability in Figure 3. The agreement is quite 
good in view of the fact that the typical distances used on the lattice remain modest. 

The same comparison has been made for the converse situation, in which the real axis is closed, except 
for the interval [21,^2] which is open. The corresponding probability 

Pr'(zi, Z2; 2) -Pi = ^ (/i'^^''(^i) ^''''^(^2) 0(2)>UHP , (40) 
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may be obtained from the previous one by the inversion z — > — l/z, which exchanges the boundary condition 
around with that around infinity. The function is the ratio of a 4-point and a 2-point function, which 
are both invariant under a global conformal transformation, so that the inversion itself has no effect at all. 
However, the function has a non-trivial monodromy around x = 1, so that the inversion interchanges the 
two behaviours H38(l around R — Q and R = +cx), implying that it must include a change of sign of all factors 
^/l — X. One finds that 

Pr'(-i^,i?;.)-A = A^df^ (41) 

is the opposite of P^^ {—R, R; z). We observe that this probability is given by the correlator (|40|l through 
the uj channel in the fusion of the two /i's, since the field does not couple to the identity. As noted before, 
this is a remnant of the regularization used to compute probabilities when the boundary is closed. 

A numerical calculation of this function on the lattice has been carried out for the same range of 
parameters as in the previous case. The results were fully consistent with the predicted change of sign. 

5. Isolated dissipation 

In this section and the next one, we introduce dissipation (or mass) at an isolated site located in the bulk 
of the UHP or on a closed boundary. The amount of dissipation does not really matter, so we consider a 
minimally dissipative site, with unit mass. It is convenient to start with the case where a single boundary 
condition is imposed along the boundary. The case where two different boundary conditions coexist is treated 
in the next section. 

^From our discusion in Section 2, introducing dissipation at a given site, located at z say, corresponds to 
raising the toppling threshold of that site from 4 to 5. This has the consequences that the height at z takes its 
values in {1,2,3, 4, 5}, and that each time it topples, one grain of sand is dissipated. In terms of the toppling 
matrix, the introduction of dissipation at n sites Zk corresponds to going from A to A + Dz^ + . . . + Dz,^ , 
where Dz = Si^z 5j^z is the defect matrix which makes the site z dissipative (one could say open). 

One can again measure the effect of introducing dissipation by computing the fraction by which the 
number of recurrent configurations increases, given by the ratio of partition functions, 

F{zu ■■■,zn)= '^''^^^''detA ""'''"^ = + ^\,ei......,z,.,. (42) 

and which reduces to the calculation of an n-by-n determinant. 

This fraction F{zi, . . . , z„) is related to the probability that all sites of the lattice have a height smaller 
or equal to 4, or to the probability that one of the sites Zi has a height equal to 5, 

F{zi, z„) = p^^^^j^^ h,<4]^l~ Proh[hz, = 5 V . . . V /i^„ = 5] ' ^^^^ 

For n — I, and when the boundary condition on the real axis is all open, the value of F{z) can easily 
be worked out for z far from the boundary, 

Fop{z) = 1 + A-pi(z; z) = ^ log |z - zH- ^ (7 + ^ log2) + 1 + . . . (44) 

where the ellipses stand for corrections that go to when \z — z\ goes to infinity (correction terms to scaling), 
and 7 = 0.57721 ... is the Euler constant. 

For the all closed boundary condition, the fraction Fc\{z) diverges, as do all higher point functions 
Fci{zi, . . . , Zn). Indeed since the sites Zj are the only sink sites, the relaxation process will produce a 
constant fiow of sand towards them, making at least one of them almost always 'full', that is, Proh[hzi = 
5V. . ■\^hz^ = 5] = 1. If one regularizes the fraction Hke in Section 3, by picking the coefficient of 2A~j^jjg(0; 0) 
in Fci{z) = 1 + A^^(z; z), then the regularized fraction is equal to 1. 
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These two simple calculations are consistent with identifying the introduction of dissipation at z with 
the insertion of a logarithmic field oj{z^z). As recalled in Section 2, this is also supported by the fact that 
J d^2;a;(z, z) is the perturbation term that drives the conformal action to a massive regime, when dissipation 
is added at all sites of the lattice. Then, 

(w(z,z))op = -!-log|z- z| +70, (w(z,z))ci = 1, (45) 

ZTT 

with7o = ^(7+|log2) + l. 

For n = 2, the calculation of F{zi, Z2) for an all open boundary yields, in the scaling limit, what should 
be the 2-point function 



{uj{zi)uj{z2))op 



1 

47r 



^log- 



2 Fi 



Z2\ 



\Zl - Z2\ 



+ ( log \zi - zil + 70 ) f 7^ log \z2 - ^2! + 70 ) • (46) 



2tt 



\2tt 



When the boundary is all closed, the regularized fraction gives 

{Uj{zi)uj{z2))cl = -log|zi2| + 270 + -5-log- 

TT ZTT \Zl 



\Zl - Z2\ 



(47) 



Zl\\Z2 - Z2\ 

These correlators allow to compute the bulk OPE of w with itself. Assuming that it transforms like 

dw 
dz 



uj{z,z) uj(w,w) + Alogl^p (its normalization is fixed by H45|l '). then the Mobius invariance fixes the 
general form of the OPE 



Lli{zi)uj{z2) 



2Alog|zi2p uj{z2)+ 5 + aAlog|zi2p-A2log2|zi2p 1(^2) + 



(48) 



for two constants a, h and where the dots stand for terms which vanish when zi — Z2. A simple comparison 
with H4()|) then yields, using H45|l . 

A-—, a = 270, & = -— = -72. (49) 



1 , a- 

Then the expectation value of the OPE in front of a closed boundary exactly reproduces 147|) . up to a 
term which vanishes if zi — Z2. This not only provides a consistent check on the OPE, but also of the 
regularization prescription we have used throughout for an all closed boundary. 

The problem of identifying the field corresponding to the insertion of dissipation at a boundary site, and 
the corresponding boundary OPE, may be analyzed along the same lines. We consider a closed boundary 
only, as sites on an open boundary are already dissipative. 

That field must belong to the representation TZq of the c — ~2 theory, and anticipating a little bit, it is 
not difficult to see that it is a weight zero logarithmic field, which we call ujb (hence a logarithmic partner 
of the identity). Indeed the regularized fractions F(xi, . . . ,a;„), for all Xi on the boundary, are supposed to 
converge in the scaling limit to the n-point function {uJb(xi) . . . W{,(a;„))ci, and one finds, for n = 1, 2, 3, 

{uJb{x))cl = 1, 

2 , , , . 5 



{uJb{xi)uJb{x2))c\ 
{LL>bixi)uJb{x2)(^b{x3))cl 



- log|a:i2| 

TT 



-470-3, 



(50) 
(51) 



log \x 



121 



270- 



1 , |2;i3 

- log 

TT X23 



/ 2 5\ / 2 

l^- log 1X12 I +470 - 2) (- l0g|xi3| 



470-2 



(52) 



They univoquely fix the first terms in the OPE of uJb with itself, given on general grounds by a chiral version 
of (EHl), 

LLib{xi)L0b{x2) ^ ab-2A;,log|a;i2| ^6(^2)+ fob+abA;, log |xi2|-A^ log^ |xi2| lAr(a;2) + . . . (53) 
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A straightforward comparison with the 1- 
coefRcients 



2- and 3-point functions yields the value of the three 



Ab — , 

TT 



h - 



(54) 



The equahty Afc = \n (see Eq. (|2fi|l l is not a coincidence. The cylinder Hilbert space with closed 
boundary condition on both edges contains a single copy of the representation TZo jnj- Therefore the two 
fields LOb and ojn must be (almost) proportional, and since their normalizations are identical, their conformal 
transformations must be the same. In fact, we will see in the next section that they are not quite identical, 
but rather differ by a multiple of the identity In, which explains why we gave them different names. On the 
other hand, there is only one possible field in TZq which transforms like the identity I^r, so that there is no 
ambiguity for the identity term. 

Finally one may examine how the bulk dissipation field uj{z,z) close to a boundary expands on boundary 
fields. It must expand on fields of Vo close to an open boundary and on fields of TZq close to a closed boundary. 
Moebius invariance again fixes the precise form of the first terms in the expansion, for z — x + iy, 



uj{z,z) — cujb{x) + + (cAfc — 2A) log2y I{x) + 



(55) 



47r 



and Xh — 



The identity I(x) can either be l£){x)) or lAr(x) depending on the type of 



where A = 
boundary. 

If it is an open boundary, the coefficient c is equal to 0, and the 1-point function (|45|) readily gives 
d — ^Q. In particular, one sees that the field corresponding to the addition of dissipation at an open site must 
be a descendant of the identity . It has been identified in as the primary weight 2 field proportional 
to dOdO in the lagrangian realization Q). 

If u;(z, z) is close to a closed boundary, the simplest way to determine the coefficients is from the 2-point 
function l|47() in the limit where zi and Z2 approach the boundary. One easily finds c = 1 and d = f — 7o- 



Summarizing, one has 

u){z, z) = Wb(x) 



r5 

4 



70 ■ 



277 



log 2y 



In{x) 



70 



27r 



log 2y 



Id{x) 



(56) 



6. Dissipation v^rith change of boundary condition 

In order to probe finer effects of the insertion of isolated dissipation, we consider, as in Section 3 and 4, the 
cases of an open boundary with a closed interval / = [zi, Z2], and the inverse situation, and the corresponding 
fractions i^op(^i, Z2; z) and Fci(2:i, Z2\z). They measure the effect of inserting a unit of dissipation at site z 
in presence of two boundary conditions, and are still given by an expression like (|42|l . where A = Aop — Bi 
or A = Aci + Bi, 

p , , det[I-A-pi(B,-i?,)] 

Fop{zi, Z2; z) = 



Br 



Fci{zi,Z2;z) 



det I 



det I 



(57) 



det [I - Aop 

The numerical calculation of these two functions is fairly straightforward. For an interval / = [—R, R] 
symmetric around the origin, we have computed the values of _Fop(— i?, i?; z) and Fci{—R, i?; z) (regularized 
as before) for z on vertical and horizontal slices in the UHP, and for different values of R (same as in Section 
4). The results are given in Figure 4 and 5, as colour dots. In both figures, the plots on the left correspond 
to the vertical slice ^z — and for three values R — 10, 30, 50, while the plots on the right correspond to 
i? = 30 for the three horizontal slices 3z = 15, 40 and 70. 

In the conformal theory, the fractions F{zi, Z2; z) should correspond to a 3~point function in the UHP. 
The two jumps of boundary condition on the real axis are effected by the insertion of two h — —■^ primary 
fields, namely one ^^'^ and one /i^'^, whereas the insertion of dissipation is represented by the insertion 
of the logarithmic field w(z, z) (see previous section). 
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Figure 4. Values of Fop{—R,R;x + iy) as function of y (left) or as function of x (right), either 
computed numerically (dots) or from CFT computations, Eq. II64II (solid lines). The curves on the 
left correspond to a; = and 5 < y < 200, for values of R equal to i? = 10 (light blue), 7? = 30 
(blue) and R = 50 (violet). Those on the right correspond to 7i = 30 and y — 15 (light orange), 
y = 40 (orange) and y = 70 (red), with -100 < x < 100. 



2 . 8 




Figure 5. Values of Fci{—R, R;x + iy) as function of y (left) or as function of x (right), numerically 
(colour dots) or from CFT calculations (solid lines). The curves on the left correspond to a; = 
and 5 < y < 200, for values of R equal to 7? = 10 (light blue), i? = 30 (blue) and R = 50 (violet). 
Those on the right correspond to i? = 30 and y — 15 (light orange), y = 40 (orange) and y = 70 
(red), with -100 < a; < 100. 



Then the scaling limits of the lattice fractions are given by 3-point correlators on the UHP or 4-point 
chiral correlators on the plane 

F{zi,Z2; z) = ^ {^J'{Zl) fJ.{z2) t^{z) t^(2;*))chiral,planc , (58) 

where the fi fields are chosen according to the case we want to consider. We first compute the general form 
of the 4-point function, and then choose the particular solution which suits the boundary conditions we 
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impose on the real axis. 

Because uj is the logarithmic partner of a primary field (the identity), the 4-point function requires a 
multiple step calculation. Indeed it satisfies an inhomogeneous differential equation (see for instance 
where the inhomogeneity depends on the correlators where each logarithmic field u is in turn replaced by 
its primary partner AI, with A = — One finds from the 3-point function 

(/x(zi)^(z2)tj(z3)) = P + a\ log , (59) 

the general form of the chiral 4-point correlator, with the same cross ratio x = ^i^lM. as before. 



(Ai(zi) ^l{z2) w(z3) uj{zi)) = zl!^^ {l~5 log \ — - - aX^ ( log | , — 

2 

+ (3X log 2/3A log ^34 + aX^ log log ] , (60) 

where a, f3, 7, S are arbitrary constants. 

For the case "closed interval in an open boundary" , the function 

Fop(zi, Z2; z) = -i^ (A*^^^(zi) A*^'^(Z2) ^(z) ^(z*)) , (61) 

can be determined by choosing the solution H60|l which has no logarithmic singularity when zi,Z2 go to 0, 
and which reproduces the fraction Fop{z) given in (|44|) . 

The first condition forces a — and 6 — 2(3X (remember that ^/l — x goes to +1), whereas the second 
one imposes 

P = A, 7 = A(7o + -log2). (62) 

TT 

This provides a completely explicit expression for Fap{zi, Z2', z), 

- log 2 log ^ H 

TT 27r y'l — x 27r 

For zi = —R and Z2 = R, and for z = x + iy, it reduces to 

Fop(-i?, i?; z) ^ 1 + -1 (7 + I log2) - log ( ^^^1 . (64) 

27r 2 ZTT I y ) 

For the cases worked out numerically and discussed above, the previous formula yields the solid curves 
shown in Figure 4. 

The expansion of Fop{—R, R; z) in the two regimes R small and R large read 

f 1 + (7 + ^ log2) + logy + . . . for i? « 1, 

F,p{^R,R;z)=\ 27r 2 2n ^ ^g^^ 

1 + (7 + o log2 + logi?2) - — logy + . . . for R » 1. 

The large R limit corresponds to shrinking the open portion to nothing, and therefore to fusing the two [i 
fields at infinity. From Hll|l . the fusion gives rise to two channels, proportional to logi? and to 1. 
The fraction for the opposite case, that of a closed boundary containing an open interval, 

i^el(^l, Z) = (m'^'^(zi) li'''''{z2) 0J{Z) iv{z*)) , (66) 

Azi2 



F,p{z^,z2;z) = -f^ + -\og2-—\og^ ^ _ ' + — log|z-z*|. (63) 
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is related to the previous one by the inversion z —1/z. Since the correlators are computed so as to be 
invariant under the (inhomogeneous) Moebius transformations, the only change comes from the monodromy 
properties discussed in Section 4, which simply change the sign of ~ x. One thus obtains 



1 . , 1 



1 



Fei(zi,Z2;z)=7o + -log2-— log ^ ' + log|z-z*|. (67) 

TT ZTT y 1 — X 2tT 

Setting zi = — i? and Z2 = R, we obtain the formula, 

1 7 1 r 1 + im-z^\ ^ 

Fei(-i?, i?; z) = 1 + — (^ + - log2) - — log { ^-^1 , (68) 

which has been used to generate the solid lines pictured in Figure 5. The agreement is excellent. The 
expansion of Fc\{~R, R; z) in the two extreme regimes now read 

l + ^(7 + ^og2-logi?2)-i^log^ + ... fori?«l, 
Fci(-i?,i?;z)=<; ^ ^ (69) 

l + — {l + -\og2) + —\ogy + ... iovR^l. 

We finish this section with two comments, and first on the i? ^ 1 limit in the previous expression. 

According to the fusion of two /i fields on a closed boundary, the logarithmic term — ^ log R should be 
proportional to Xn{uj{z, z))c\ logi? and confirms all previous results. On the other hand, the non-logarithmic 
piece should correspond to (wAr(O) lo{z, z))ci, namely a chiral 3-point function on the plane (a;jv(0) lli{z)lo{z)) . 
Its z dependence is however unusual and is due to the fact that the two logarithmic fields involved have 
different inhomogeneous terms in their conformal transformations, Aat 7^ A. An explicit calculation shows 
that the 3-point function in this situation is exactly what the above limit yields. 

Our second remark is on the relation between the two logarithmic boundary fields, ujn and ujb- Both have 
the same normalization {u!t,{x)) = {u!n{x)) = 1, and the same conformal transformations Xn = Af, = — — , 
yet they do not have the same 2-point function. From the limit 2:12,2:34 ^ of the four jj correlator H27|l. 
one obtains 

2 4 

{lun{xi)lun{x2))c\ = -log|a::i2| + -log2, (70) 

TT TT 

which differs from {LUi,{xi)u!b{x2)) in (|51|l by the constant piece. The two fields belonging to the same 
representation, they can only differ by a multiple of the identity. Comparing the two 2-point functions, one 
finds 

3 1 1 

LJb^LJN + KlN, K = - + -(7 - -l0g2). (71) 

4 TT Z 

One obtains the same value of k from the mixed 2-point function arising in the limit Z12 0, z — z* ^ 
of the correlator ()66|l . when using the field decomposition (|56|) . 



7. Conclusions 



Our primary motivation in this article was to contribute to establishing the Abelian sandpile model as 
a model that can be described by a conformal field theory. Whether and to what extent the conformal 
description holds remains an unclear issue, mainly because of the highly non-local interactions present in 
the sandpile model. 

We have essentially considered three operators. The first is the operator that changes a boundary 
condition between open and closed; the second corresponds to the height 1 variable; and the third one is 
the insertion of dissipation at a non-dissipative site. We have studied fine details of their mixed correlations 
on the upper half plane, and found indeed a full agreement with the predictions of a logarithmic conformal 
theory with central charge c — —2. 
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Though our results are certainly encouraging and show that the conformal description cannot be all 
wrong, it is still however far from proving that every aspect of the sandpile model has a counterpart in the 
conformal theory and vice-versa. Before this goal may be attained, important questions must be answered, 
like the existence of other boundary conditions than open and closed, the sandpile interpretation of the 
weight 3/8 primary field, the sandpile significance of the Z2 symmetry and of the W-algebra present in the 
conformal theory, the conformal description of the higher height variables and of the avalanche observables. 
The list is not exhaustive but only shows the benefit one can expect on both sides. 



Appendix 

We give here some details on the numerical calculations reported in the text. All numerical results are 
related to the calculation of finite determinants, the largest ones being of a typical size of 100, though larger 
ones have been considered. 

The matrices of which the determinant is to be computed are of the form (I+A^p'^A) or (I+A^j^A) where 
A is a numerical matrix. Their entries are linear combinations of (Aop)jj^ and (Ad)^^, where i, j e Z x Z> 
run over some set of lattice sites in the upper half-plane, the boundary of which is the horizontal line y = 1. 
By the method of images, the inverse matrices A^p'^ and A^^ are related to the inverse Laplacian on the full 
plane Apj^^^^^. If (toi,71i) and (7712,712) are the integer coordinates of i and j respectively, then 

^^°p)(mi,ni),(mi,n'2) ^ (^P^'^'^°\mi,ni),{m2,n2) ~ (^P^'^'^<'\mi ,ni) ,{m2 ,-112) ' (^-1) 
(^cl)(mi,„i)^(m2,n2) ~ (^Planc) (mi ,ni ) ,(™2 ,"2) '■^P''^"<= ) ("H ,rii),()n2 ,l-n2) ' (^'2) 

Thus the entries of the inverse Laplacian on the plane 7? are required. Because of the horizontal 
and vertical symmetries, it is enough to know (Apianc)ij^ where one site is the origin. A simple Fourier 
transformation yields a divergent integral representation 

ff^^ d^k Qikim+ik2n 
(Aplanc)^i,„)^(0,o) = JJ^ ^ 4-2c0sfci-2c0Sfc2' ^^'^^ 

and a finite representation for the difference 

ff^^ d^k Qikim+ik2n _ ^ 

(/.(m, n) = (Apla„e)^i,„),(o,0) - (^Planc)ro'o),(0,0) = Jj^ ^ 4 - 2 COS fcl - 2 COS fca ' ^^'^^ 

So the entries of A^p'^ are well-defined and finite. On the other hand, the entries A^j^, being a sum of two 
inverse Laplacian entries, have all an infinite piece, which can be taken as 2(Apianc)(o^Q) (q o)' "^^^ regularized 
determinant Q is however finite as its entries are differences of (Ad)^^. 

^From the reflection symmetries, (j){m,n) = (j)(—m,n) = (f>{m,—n) = (j){n,m)^ it is enough to known 
0(771, 77) in the first half-quadrant delimited by the positive horizontal axis m > 0, 77 = 0, and the diagonal 
m — n> 0. A convenient procedure ^Bj to compute (j){m, n) in that region is to use the known exact values 
on the diagonal, given by 

1 1 

0(7^,777)---^^^—^, (771 >1) (A.5) 

k=l 

and then to propagate the function (p from the diagonal down to the x-axis by a repeated use of the Poisson 
equation, 

4(j){m, n) — (j){m + 1, 71) — (f>{m — 1, 77) — ipim, 71 + 1) — (f>{m, 77 — 1) = Sm.o Sn,o- (^-6) 

In this way the knowledge of 0(777, 77) for fixed 777 and for < n < 777, and of the diagonal element 
0(777 + 1, 777+1), allows to determine all the values of on the next line, namely 0(777+1, 77) for < 71 < 7?7 + l. 
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This way of propagating the function <j) is nuniericaUy unstable because the Poisson equation involves 
the difference of close numbers. If the propagation is not performed with enough numerical precision, the 
resulting values of (j) depart very wildly from what is expected. In the computations reported in the text, 
the values of 4'{m, n) are required for values of m of the order of 400. For the above propagating procedure 
to produce sensible results, all calculations were performed on 320 decimal places. 

Once the actual values of </)(to, n) are obtained, the determinants can be computed. All determinants 
considered in the text diverge exponentially, or go to exponentially, with their size. Consider for instance 
det[I — A~p-^i3/], where / is an interval on the boundary, possibly disconnected (like in (|^ ). In the sandpile 
model, it is equal to the number of recurrent configurations when the boundary is open except on the 
segment / which is closed, divided by the corresponding number with an all open boundary. Because a 
closed boundary site has a free energy smaller than an open boundary site, by an amount equal to IGj-K 0, 
the determinant is dominated by an exponentially small term e^^'^'^'/'^ [G is the Catalan constant). For the 
same reason, the determinant det[I + A^j^B/] for the converse situation is dominated by an exponentially 
diverging term e^'^\^\l'^ . The same is true if a matrix (relevant for a unit height) or (used for an 
isolated dissipative site) is added to B/. 

These exponential factors drop out in ratios of partition functions eventually related to 4-point CFT 
correlators — namely the ratios in Eqs. H23|l. (|28|l . (|31|l and H57|l — , but taking the ratio of huge numbers 
is not numerically efficient. To avoid this problem, one multiplies the matrices by the proper factor e^'^'^l'^ 
before computing its determinant, so as to kill the dominant exponentials. 

As the determinant calculations generate a moderate loss of precision, the precision on the matrix entries 
is at this stage lowered to 25 decimal places. The numerical errors on the final results are expected to be 
smaller than 0.001%. In view of the relative importance of the corrections to scaling, there is no need to 
improve it. 
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